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We demonstrated microwave-induced cooling in a superconducting flux qubit. 
The thermal population in the first-excited state of the qubit is driven to a 
higher-excited state by way of a sideband transition. Subsequent relaxation 
into the ground state results in cooling. Effective temperatures as low as Tcs ~ 
3 millikelvin are achieved for bath temperatures Tbath = 30 — 400 millikelvin, 
a cooling factor between 10 and 100. This demonstration provides an ana- 
log to optical cooling of trapped ions and atoms and is generalizable to other 
solid-state quantum systems. Active cooling of qubits, applied to quantum in- 
formation science, provides a means for qubit-state preparation with improved 
fidelity and for suppressing decoherence in multi- qubit systems. 
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Cooling dramatically affects the quantum dynamics of a system, suppressing dephasing and 
noise processes and revealing an array of lower-energy quantum-coherent phenomena, such as 
superfluidity, superconductivity, and the Josephson effect. Conventionally, the entire system un- 
der study is cooled with ^He-^He cryogenic techniques. Although this straightforward approach 
has advantages, such as cooling ancillary electronics and providing thermal stability, it also has 
drawbacks. In particular, limited cooling efficiency and poor heat conduction at millikelvin 
temperatures limit the lowest temperatures attainable. 

A fundamentally different approach to cooling has been developed and implemented in 
quantum optics (17]|2]|5]|?1). The key idea is that the degrees of freedom of interest may be cooled 
individually, without relying on heat transfer among different parts of the system. By such di- 
rected cooling processes, the temperature of individual quantum states can be reduced by many 
orders of magnitude with little effect on the temperature of surrounding degrees of freedom. In 
one successful approach, called sideband cooling (l5]|^0|8|), the unwanted thermal population 
of an excited state 1 1) is eliminated by driving a resonant sideband transition to a higher excited 
state 1 2), whose population quickly relaxes into the ground state |0) (Fig. [IK)- The two-level 
subsystem of interest, here {|0), |1)}, is efficiently cooled if the driving-induced population 
transfer to state |0) is faster than the thermal repopulation of state The sideband method, 
originally used to cool vibrational states of trapped ions and atoms, allows several interesting 
extensions ([71l2ll3]|?ll9l [7Ql[77|[72]) . For example, the transition to an excited state can be achieved 
by nonresonant processes, such as adiabatic passage @, or by adiabatic evolution in an optical 
potential <\10\\ll\U2\i . Other approaches, such as optical molasses and evaporative cooling, have 
been developed to cool the translational degrees of freedom of atoms to nanokelvin tempera- 
tures, establishing the basis for the modem physics of cold atoms (1731) . 

Superconducting qubits are mesoscopic artificial atoms (f7?|) which exhibit quantum-coherent 
dynamics (I75l) and host a number of phenomena known to atomic physics and quantum op- 
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tics, including coherent quantum superpositions of distinct macroscopic states (|76][77l) . time- 
dependent Rabi oscillations (I7^[79ll20]l271l22l[23ll24)) . coherent coupling to microwave cavity 
photons (l25l 12611271) and Stiickelberg oscillations via Mach-Zehnder interferometry (|28l 12911301) . 
In a number of these experiments, qubit state preparation by a dc pulse or by thermalization 
with the bath was used. It is tempting, however, to extend the ideas and benefits of optical 
cooling to solid-state qubits, because they present a high degree of quantum coherence, a rel- 
atively strong coupling to external fields, and tunability, a combination rarely found in other 
fundamental quantum systems. 

We demonstrate a solid-state analog to optical cooling utilizing a niobium persistent-current 
qubit (|371) . a superconducting loop interrupted by three Josephson junctions (l32l) . When the 
qubit loop is threaded with a dc magnetic flux /q ^ $o/2, where $o = h/2e is the flux 
quantum Qi is Plancks constant), the qubit's potential energy exhibits a double-well profile 
(Fig. [IK), which can be tilted by adjusting the flux detuning, 5/q = /q — $0/2, away from 
zero. The lowest-energy states of each well are the diabatic qubit states of interest, |0) and 
characterized by persistent currents Jq with opposing circulation, whereas the higher-excited 
states in each well, e.g., |2), are ancillary levels that form the "sideband transition" with the 
qubit. In contrast to conventional sideband cooling, which aims to cool an "external" harmonic 
oscillator (e.g., ion trap potential) with an "internal" qubit (e.g., two-level system in an ion), our 
demonstration aims to cool an "internal" qubit by using an ancillary "internal" oscillator-like 
state [supporting online material (SOM) Text]. 

When the qubit is in equilibrium with its environment, some population is thermally excited 
from the ground-state |0) to state |1) according to pi/po = exp[— (^i — £:o)/^B^bath], where 
Po,i are the qubit populations for energy levels £0,1' is the Boltzmann constant, and Tbath is 
the bath temperature. To cool the qubit subsystem below Tbath, in analogy to optical pumping 
and sideband cooling, a microwave magnetic flux of amplitude A and frequency v targets the 
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Figure 1 : Sideband cooling in a flux qubit. (A) External excitation transfers the thermal population from state 
|1) to state |2) (straight line) from which it decays into the ground state |0). Wavy lines represent spontaneous 
relaxation and absorption, ^ r2i, Fqi. The double well is the flux-qubit potential comprising energy levels. 
(B) Schematic band diagram illustrating the resonant and adiabatic sideband cooling of the ac-driven qubit. 1 1) ^ 
1 2) transitions are resonant at high driving frequency i/ (blue lines) and occur via adiabatic passage at low u (red 
lines). Agi and A12 are the tunnel splittings between |0) - |1) and |1) - |2). (C) Cooling induced by ac -pulses 
with driving frequencies ly = 800, 400, 200 and 5 MHz. State |0) population versus flux detuning <5/q for the 
cooled qubit and for the qubit in thermal equilibrium with the bath (black lines, Tbath = 300 mK). Measurements 
for 1/ — 800, 200 and 5 MHz are displaced vertically for clarity. (Inset) Pgw versus Sfq over a wider range of flux 
detuning; j/ = 800 MHz. 

|1) 1 2) transition, driving the state- 11) thermal population to state |2) from which it quickly 
relaxes to the ground state |0). The hierarchy of relaxation and absorption rates required for 
efficient cooling, ^ , Fqi, is achieved in our system owing to a relatively weak tunneling 
between wells, which inhibits the inter-well relaxation and absorption processes, |2) 1 1) and 
|0) compared with the relatively strong intra- well relaxation process |2) |0). This 

three-level system behavior is markedly different from the population saturation observed in 
two-level systems. 

The cooling procedure illustrated in Fig. [TK is generalized to the energy-band diagram 
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shown schematically in Fig. [T^. The diabatic-state energies, 

£1,0 = ±V/q, e2',2 = e* ± Ic^Sfc^, (1) 

are linear in the flux detuning 5/q, with the energy e* ~ 25 GHz and Jq = 1.44 GHz/m$o in 
our device, and exhibit avoided crossings Aqi ~ 12 MHz and A12 = A02' ~ 100 MHz due 
to quantum tunneling through the double-well barrier (Fig. [IK)- The diabatic levels exchange 
roles at each avoided crossing, and the energy band is symmetric about 5/q = (l53l) . 

Under equilibrium conditions, the average level populations exhibit a thermally-broadened 
"qubit step" about 5fc^ = , the location of the |0) - |1) avoided crossing. This is determined 
from the switching probability Pgw of the measurement superconducting quantum interference 
device (SQUID) magnetometer, which follows the |0) state population (1321) . 

1 £ 

Psw =0(1 + ^"^0), mo = tanh 77^-^, (2) 

where F is the fidelity of the measurement, mo = Po — Pi is the equilibrium magnetization that 
results from the qubit populations po,i, T = Tbath, and £ = ei — £0 oc 5fq as inferred from 
Eq. [TJ In the presence of microwave excitation targeting the |1) |2) transition, the resultant 
cooling, which we will later quantify in terms of an effective temperature T^s < Tbath^ acts to 
increase the ground-state population and, thereby, sharpen the qubit step. This cooling signature 
is evident in Fig.[l]I!, where we show the qubit step before and after applying a cooling pulse at 
several frequencies for Tbath = 300 mK. 

The cooling presented in Fig. [H B and C, exhibits a rich structure as a function of driving 
frequency and detuning, resulting from the manner in which state |2) is accessed. The |1) — > 
1 2) transition rate can be described by a product of a resonant factor and an oscillatory Airy 
factor (l30l) . The former dominates at high frequencies (800 and 400 MHz), where well-resolved 
resonances of n-photon transitions are observed, as illustrated in Fig. [iB (transition in blue) and 
Fig.[T]C (top traces and inset). The cooling is thus maximized near the detuning values matching 
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e2—S\ = nhu (downward arrows in Fig.[T]Z!). At intermediate frequencies (400 and 200 MHz), 
the Airy factor becomes more prominent and accounts for the Stiickelberg-like oscillations that 
modulate the intensity of the n-photon resonances (|2^|30|) . Below v = 200 MHz, although 
individual resonances are no longer discernible, the modulation envelope persists due to the 
coherence of the Landau-Zener dynamics at the A12 avoided crossing (|3Q|) . The |1) |2) 
transition becomes weak near the zeros of the modulation envelope, where we observe less 
efficient cooling, or even slight heating (e.g., upward arrows in Fig. \Vp, 800 and 400 MHz). 
This is a result of the |0) — > |1) transition rate which, although relatively small, Aq^^ <^ AI2, 
acts to excite the qubit when the usually dominant |1) |2) transition rate vanishes. At low 
frequencies [z^ < z^o = (^12/^*)^^^ ~ 10 MHz], the state |2) is reached via adiabatic passage 
(Fig. [T^, red lines) and the population transfer and cooling become conveniently independent 
of detuning (see u = 5 MHz in Fig. [T]C). 

Maximal cooling occurs near an optimal driving amplitude (Fig. |2l). Fig. [2K shows the |0) 
state population Psw measured as a function of the microwave amplitude A and flux detuning 
6fq for frequency u = 5 MHz. The adiabatic passage regime, realized at this frequency, is 
particularly simple to interpret, although higher frequencies allow for an analogous interpreta- 
tion. Cooling and the diamond feature of size A* = £* /2Iq can be understood in terms of the 
energy band diagram (Fig. [IB). For amplitudes < A < A* /2, population transfer between 
states |0) and |1) occurs when A > \5fq\, such that the sinusoidal flux reaches the Aqi avoided 
crossing; this defines the front side of the observed spectroscopy diamond symmetric about the 
qubit step. For amplitudes A* /2 < A < A*, the A12 (A02') avoided crossing dominates the 
dynamics, resulting in a second pair of thresholds A = A* — \5fq\, which define the back side 
of the diamond. As the diamond narrows to the point A = A*, cooling is observed. There only 
one of the two side avoided crossings is reached and, thereby, strong transitions with relaxation 
to the ground state result for all 5fq, yielding the sharpest qubit step. For A > A*, both side 
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Figure 2: Optimal cooling parameters and effective temperature. (A) State |0) population Pgw versus flux 
detuning, -A* /2 < dfq < A*/2, and driving amplitude A with = 5 MHz, tp = 3^s, and Tbath = 150 mK. 
Optimal conditions for cooling are realized at A = A*, where A* is defined in Fig. IB. (B) Effective temperature 
Tcff. Qubit steps measured at = 5 and 245 MHz (circles) and best fits to Eq. |2] At 245 MHz, the aggregate 
temperature fitting (blue, Tcff = 13.8 mK) and the convex fitting (orange. Toff = 9.1 mK) are shown. Tbath = 30 
mK. (C and D) (Inset) Detail of the region A A* [white box in (A)] for = 5 MHz (top) and = 245 MHz 
(bottom). In each case, T^ff is extracted from the qubit step as in (B). Lines are guides for the eye; tp = 3 /is, 
Tbath = 150 mK. (E) Spectroscopy of uncooled (top) and cooled (middle) qubit (5 MHz, 3-/is cooling pulse) at 
Tbath = 30 mK. Cumulative switching -probability distribution as a function of /g and Sfq under 500-MHz ac 
excitation. 

avoided crossings A12 and A02' are reached simultaneously when \5fq\ < A — A*, leading again 
to a large population transfer between |0) and 1 1). 

When an ac field is applied, the qubit is no longer in equilibrium with the bath, but it can 
still be well-characterized by an effective temperature Teg using Eq. [21 with T = Teg. This is 
illustrated in Fig. |2^ for z/ = 5 MHz and u = 245 MHz (Tbath = 30 mK). At z/ = 5 MHz, 
the qubit step clearly follows Eq. [2l as shown with a fit line for Teg = 3.4 mK. At 245 MHz, 
individual multiphoton resonances are evident, and Psw is a non-monotonic function of 5fq. 
In this case, T^s is still a useful parameter to quantify the effective cooling, but it should be 
interpreted as an aggregate temperature over all frustrations. Alternatively, because the cooling 
is maximized at individual resonances, one may perform a convex fitting of Eq. |2l where only 
the solid (orange) symbols are taken into account to determine the effective temperature at the 
resonance detunings. The convex effective temperature = 9.1 mK is smaller than the 
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aggregate value Tes = 13.8 mK. In the remainder of the paper, we refer to the more conservative 
effective temperature obtained using the aggregate definition. 

Figure [2l C and D, show the variation of Tgfj about A = A* for u = 5 MHz and u = 245 
MHz, respectively, in the region marked with a white rectangle in Fig. [JK (insets show the raw 
data). As seen in these figures, T^s typically presents a minimum, where the cooling is most 
efficient and from which A* can be determined. 

To determine whether the observation of a sharp qubit step proves that the system makes 
transitions to the ground state, as opposed to selectively populating an excited state with the 
same magnetization, we measured the excitation spectra of the "pre-cooled" qubit and of the 
qubit in thermal equilibrium with the bath (Fig. [2E)- In the former, a weak ac excitation was 
applied immediately after the cooling pulse (time-lag less than 100 ns), well before the system 
equilibrates by warming up to the bath temperature (see below). By comparing the excitation 
spectra of the equilibrium and cooled systems (Fig.[2E^, Tbath = 30 mK), we note that, although 
cooling markedly reduces the step width, making the qubit much colder, the excitation spectrum 
remains unchanged. Because the ac excitation is resonant with the |0) |1) transition only, 
this strongly indicates that the population in a cooled qubit is in the ground state. 

Figure m A and B, summarize the dependence of T*g = Teg{A*) on the dilution refrigerator 
temperature Tbath = 30 — 400 mK for several frequencies u, spanning the resonant sideband 
to the adiabatic passage limits, with a fixed pulse width tp = 3 yus. In Fig. [3K, at large u, T*^^ 
exhibits a monotonic increase with Tbath^ which becomes less pronounced as u decreases. In 
the adiabatic passage limit, e.g. u = 5 MHz, T*g ^ 3 mK is practically constant and reaches 
values that, notably, can be more than two orders of magnitude smaller than Tbath- In Fig - 
T*g is observed to increase linearly with u for different values of Tbath- Because the number of 
resonances in the qubit-step region is inversely proportional to u, the cooling at the individual 
resonances depends only weakly on u when using the convex definition T^^^'"^'^{A*). 
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Figure 3: Effective temperature T*g for A = A* and measurement fidelity F. (A) T*g versus Tbath at the 
indicated driving frequencies ly. T*^^ increases with Tbath at high but remains constant at low i^. (B) T*g versus 
v for different Tbath- Lines are linear fits. (C) F versus Tbath at the indicated ly. A pulse width ip = 3 /is was 
used in all cases. 

Figure [3]C displays the measurement- fidelity F versus Tbath- Although the qubit is effec- 
tively cooled, T*ff <^ Tbathj over the range of Tbath in Fig- HI A and B, the readout SQUID is 
not actively cooled, and its switching current distribution broadens with Tbath (fig- S2). At high 
temperatures, the fidelity F, defined in Eq. |2l becomes too small to discriminate the two qubit 
states; this is independent of the qubit's effective temperature, which remains ~3 mK at all 
values of Tbath- We observe that the fidelity F is larger than 0.8 for Tbath < 100 mK, remains 
above 0.5 at ^He refrigerator temperatures, but drops to T ^ 0.1 at Tbath = 400 mK, limiting 
our ability to measure the qubit state at higher temperatures (SOM Text). 

The cooling and equilibration dynamics of the qubit are summarized in Fig. H] (Tbath = 150 
mK). Cooling a qubit in equilibrium with the bath requires a characteristic cooling time. In turn, 
a cooled qubit is effectively colder than its environment, a non-equilibrium condition, which 
over a characteristic equilibration time will thermalize to the environmental bath temperature. 
This relation between cooling and equilibration times determines the facility of cooling the 
qubit and performing operations while still cold. Fig. |4l A and B, show the time evolution at 
cooling and warming up of the qubit step. The top panels show Psw as a function of 5fq and 
cooling-pulse length tp (Fig. UK, u = 245 MHz), and as a function of 6fq and waiting-time 
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Figure 4: Dynamics of cooling and equilibration. (A) (Upper panel) State |0) population P^w as a function of 
6fq and cooling pulse width tp {v = 245 MHz). (Lower panel) Teff versus tp extracted from upper panel (circles) 
and exponential fit (blue line) with ^ l-/^s time constant. (B) (Upper panel) State |0) population Psv, as a function 
of (5/q and waiting time after the cooling pulse (<p = 3 /is and v = 5 MHz). (Lower panel) Toff versus 
extracted from upper panel (circles) and exponential fit (blue line) with ^ 100-/xs time constant. Tbath = 150 
mK. (C) Characteristic equilibration and cooUng times for different Tbath- Cooling is performed at the indicated 
frequencies. 



tw after pre-cooling with a 5 MHz pulse (Fig.S^) (for tp and definition, see fig. SI). Note 
the difference in the time scales, where it is observed that substantial cooling is accomplished 
within 1 jUS (Fig. UK), but equilibration occurs over a much longer time scale (Fig.H^). Fitting 
to Eq. [2] yields Teg as a function of tp and (Fig. |4l A and B, bottom panels). The near 
exponential behavior of Teg versus tp and allows one to infer the characteristic cooling and 
equilibration times as defined by an exponential fitting (solid blue lines), which are summarized 
in Fig. lit. Notably, the cooling characteristic time is nearly independent of both u and Tbath 
and, on average, is about 500 ns. In contrast, at the base temperature of the dilution refrigerator, 
the equilibration time is about three orders of magnitude longer, 300 /iS, and remains one order 
of magnitude longer at 250 mK, a temperature that is accessible with ^He refrigerators. 
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The minimum qubit effective temperature demonstrated in this work was estimated to be 
Teff ~ 3 mK. This value is consistent with the inhomogeneously broadened linewidth observed 
in the experiment, which likely places a lower limit on the measurable minimum temperature. 
The microwave-induced cooling presented here can be applied to problems in quantum informa- 
tion science, including ancilla-qubit reset for quantum error correcting codes and quantum- state 
preparation, with implications for improved fidelity and decoherence in multi-qubit systems. 
This approach, realized in a superconducting qubit, is generalizable to other solid-state qubits 
and can be used to cool other on-chip elements, e.g. the qubit circutry or resonators 
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Materials and Methods 



Measurement Scheme The qubit consists of a superconducting loop interrupted by three 
Josephson junctions (Fig. S(5l\), one of which has a reduced cross-sectional area. DC and pulsed 
microwave (MW) magnetic fields generate a time-dependent magnetic flux f{t) = f^"^ + f^'^{t) 
threading the qubit. Transitions between the qubit states are driven by the pulsed MW flux 
/^'^(t) = A cos 2'Kut of duration tp = 10 ns— 3 /is (Fig. S(5^) with frequency v, and amplitude A, 
where A is proportional to the MW-source voltage Vac- The qubit states are read out using a DC- 
SQUID, a sensitive magnetometer that distinguishes the flux generated by the qubit persistent 
currents, Jq. After a delay t„ following the excitation, the readout is performed by driving the 
SQUID with a 20-ns "sample" current 1^ followed by a 20-;US "hold" current (Fig. S|5^). The 
SQUID will switch to its normal state voltage if h > hwfl ih > Isw,i), when the qubit is in 
State |0) (|1)). By sweeping /g and flux detuning, while monitoring the presence of a SQUID 
voltage over many trials, we generate a cumulative switching distribution function (see Fig. 2E 
and Fig. S2 below). By following a flux-dependent sample current /sw,o < h < ^sw,i we obtain 
the switching probability Ps„ that characterizes the population of state |0), and that reveals the 
"qubit step" shown in Fig. IC. 

The experiments were performed in a dilution refrigerator with a 12-mK base temperature. 
The device was magnetically shielded with 4 Cryoperm-10 cylinders and a superconducting 
enclosure. All electrical leads were attenuated and/or filtered to minimize noise. 

Device Fabrication and Parameters The device (Fig. S(5K) was fabricated at MIT Lincoln 
Laboratory on 1 50 mm wafers in a fuUy-planarized niobium trilayer process with critical current 
density Jc ~ 160 A/cm^. The qubit's characteristic Josephson and charging energies are Ej ^ 
(27r^)300 GHz and Eq ~ (27r/i)0.65 GHz respectively, the ratio of the qubit Josephson junction 
areas isa ^ 0.84, and the tunnel coupling A ^ (27r^)0.01 GHz. The qubit loop area is 16 x 16 
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fj,nr, and its self inductance IS l^n ^ 30 pH. The SQUID Josephson junctions each have critical 
current Ico ~ 2 fiA. The SQUID loop area is 20 x 20 fim^, and its self inductance is Ls ~ 30 
pH. The mutual coupling between the qubit and the SQUID is M 25 pH. 

Supporting Text 

Microwave cooling and optical resolved sideband cooling While the microwave cooling 
(MC) demonstrated in this work has similarities with optical resolved sideband cooling (RSC), 
the analogy is not a tautology. We describe here the similarities and distinctions between the 
two cooling techniques. 

In both cases, the spectra of the energy levels can be reduced to the three-level system 
defined in Fig. lA (main text). Although, as described below, the origin of the three levels is 
different for MC and RSC, cooling is similarly achieved by driving the thermal population in 
state 1 1) to an ancillary state |2), from which it quickly relaxes to state |0). 

In the case of RSC of an atom, the three-level system illustrated in Fig. lA (without the 
double-well potential) results from a two-level atomic system (TLS) combined with a simple 
harmonic oscillator (SHO) from the trap potential. Using the notation |TLS state, SHO state) 
with TLS states {\g), |e)} and SHO states {|n)} with n = 0, 1, 2, . . ., one can identify from 
Fig. 1 A the following: |0)rsc = \g,n), |1)rsc = \g,n + 1), and |2)rsc = \e,n). Note that the 
TLS represents an "internal" atomic state, whereas the SHO is an "external" trap state. 

In the case of MC of a flux qubit, the double-well potential illustrated in Fig. lA com- 
prises two coupled SHO-like wells. The left and right wells correspond to the diabatic states 
of the qubit, the clockwise and counterclockwise circulating currents, and together form the 
qubit TLS. Each well independently has a series of SHO-like states. Using the same nota- 
tion |TLS state, SHO state) with TLS states {\g), |e)}, associated respectively to the left and 
right well in Fig. lA, one can identify the following: |0)mc = IS'^O), |1)mc = |e,0), and 
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|2)mc = \9, !)• Further higher-excited states are not expUcitly shown in Fig. lA. Note that all 
states here are "internal," and that the SHO-nature of the left and right wells is limited by the 
degree of tunnel coupling between wells. 

From the above discussion, it is clear that in both RSC and MC there is a TLS combined 
with a SHO, however the roles of the levels interchange. Considering an associated frequency 
<^TLS = -E'TLs/^ for the TLS and a plasma frequency a;sHo — Esuo/K the two cases can be 
summarized as follows: 

• In the RSC case, o^sho < <^tls> and it is the SHO that is cooled. One drives population 

from the TLS ground state with higher SHO energy, |1)rsc = Ifi',^ + 1), to the TLS 
excited-state with lower SHO energy, |2)rsc = |e, n), from which it relaxes to the ground 
state |0)rsc = \g,n). 

• In the MC case, a;TLS < <^sho> and it is the TLS that is cooled. One drives population 
from the TLS excited state with low SHO energy, |1)mc = |e, 0), to the TLS ground- 
state with higher SHO energy, |2)rsc = \g, 1), from which it relaxes to the ground state 

|o)mc = b,o). 

Thus, if one considers a three-level state configuration without tagging the states with 
the terms "external" or "internal", in one cooling cycle, one cools the subsystem of interest 
{|0), |1)} by driving transitions to an ancillary state |2), which relaxes quickly to the ground 
state. For MC, the cooled subsystem is the TLS, whereas for RSC, the cooled subsystem is the 
SHO. Note that in the RSC case, because it is the SHO that is cooled, the cooling cycles can 

be cascaded to cool the multiple SHO states n = 0,1,2, However, in the MC case, it is the 

TLS that is cooled, which requires only a single cooling cycle from the TLS excited state to its 
ground state. In the MC case described in this work, all transitions are allowed, and the energy 
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levels are widely tunable. However, in many cases of RSC of atoms, certain transitions may be 
forbidden, and there is typically only limited energy-band tunability. 

Effective temperature and fidelity vs. bath temperature The qubit step and the SQUID 
switching-current distribution broaden with Tbath- The qubit can be cooled effectively, T*^^ <^ 
^bath, over the range of Tbath in Fig. 3A and Fig. 3B. However, the readout SQUID is not ac- 
tively cooled, and its switching current distribution broadens with Tbath- This is observed in 
Fig. S(6l where we plot the cumulative switching-distribution as a function of Is and 5fq of un- 
cooled and cooled qubit (5 MHz, 3-/xs cooling pulse) at different Tbath- At high temperatures, 
the switching-current distribution becomes broad and the measurement fidelity F drops to val- 
ues that become too small to discriminate the two qubit states; this is independent of the qubit's 
effective temperature, which remains about 3 mK at all Tbath - 
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Supporting Figures 




Figure 5: (Fig. SI) Schematic flux qubit and measurement scheme. (A) Qubit (inner green loop) and induc- 
tively coupled DC SQUID magnetometer (outer black loop). Josephson junctions are indicated by crosses. Red 
and blue arrows indicate the two directions of the persistent current in the qubit, /q. DC and pulsed MW magnetic 
fields f{t) = f'^'^ + f^'^{t) control the qubit and drive transitions between its quantum states. (B) The qubit state is 
inferred by sampling the SQUID with a pulsed current /g and performing voltage readout of the switching pulse, 
Vs. The current pulse Is is set to maximize the difference in the SQUID switching probability between the two 
qubit states which couple to the SQUID through the associated persistent currents, Iq. The SQUID only switches 
when the qubit is in state |0) and the switching probability Pgw measures its population. 
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Figure 6: (Fig. S2) Cumulative switching distribution of the qubit as a function of Ig and Sfq of the qubit in 
equilibrium with the bath (top) and of the cooled qubit (bottom) at different Tbath- The cooling pulse tp = 3 ^s, v 
= 5 MHz. Although the readout SQUID switching distribution broadens as Tbath increases, reducing the readout 
fidehty, the cooled qubit step remains sharp (Teff ~ 3 mK). 
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